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We investigate the sensitivity of quantum systems that 
are chaotic in a classical limit, to small perturbations of their 
equations of motion. This sensitivity, originally studied in the 
context of defining quantum chaos, is relevant to decoherence 
in situations when the environment has a chaotic classical 
counterpart. 

The exponential divergence of two trajectories, evolv- 
ing under identical equations of motion from slightly dif- 
ferent initial conditions - the famous butterfly effect - 
is a fingerprint of chaos in classical mechanics. How- 
ever, an analogous definition of "quantum chaos" based 
on evolution in Hilbert space is problematic: by unitar- 
ity, the overlap between two evolving wavefunctions - 
a natural indicator of distance between them - is pre- 
served with time, hence there is no divergence. To ad- 
dress this difficulty, Peres has suggested an alternative 
approach, in which one considers two trajectories evolv- 
ing (in phase space or Hilbert space) from identical ini- 
tial conditions but under slightly different equations of 
motion, rather than the other way around. Classically, 
even for small perturbations 0, one generically expects 
rapid divergence when the systems are chaotic according 
to the usual definition, as the perturbation (i.e. the dif- 
ference between equations of motion) soon introduces a 
small displacement between the trajectories. Quantally, 
the overlap between the wavefunctions begins at unity, 
then decays with time, and Peres suggested the rate of 
this decay - a measure of the sensitivity of quantum evo- 
lution to perturbations in the equations of motion - as a 
signature of quantum chaos. 

The sensitivity of quantum evolution also plays an im- 
portant role in the context of environment-induced de- 
coherence [HJ|]. As an illustrative example, consider a 
composite system consisting of a two-state spin (s) and 
a generic environment {£), governed by a Hamiltonian of 
the form 

H = H £ ®I S + V?®\+)(+\. + Vf<S\-)(-\ a . (1) 

Here the identity I s and the two projection operators 
act on the Hilbert space of the spin, whereas V£~ and 
Hs act on that of the environment. We view Hg as the 
"bare" Hamiltonian for the environment, and VjT as a 
perturbative coupling to the state of the spin. An initial 

state *(0) = ips(0) a\+) s + 0\—) 3 evolves into 



*(t) = e^+(t)|+).+/tyfi(t)|-). 



(2) 



pure state of the spin, a\+) s + [3\—) s , eventually becomes 
a mixture of the pointer states {|+) s , |— )s} as a result of 
monitoring by the environment. The decay of | (tpg \ipg ) | 2 
is an indicator of this process: once this overlap becomes 
negligible, the state of the spin alone can be described 
in terms of classical probabilities rather than quantum 
amplitudes. 

In view of these considerations, we are motivated to 
ask, what limits are placed on the sensitivity of a quantum 
system to perturbations in its equations of motion? The 
aim of this Letter is to provide answers to this question, 
with emphasis on systems that are chaotic in the classical 
limit. The object of our considerations will be a pair 
of wavefunctions, ipi(t) and ip2(t), identical at t = 0, 
that evolve under slightly different Hamiltonians, Hi and 
H2 = H± + V, respectively. Our measure of "sensitivity" 
will be the rate of decay of the overlap 



o q {t) = \{Mt)\Mt)W 



(3) 



We will first derive, from the uncertainty principle, a 
quite general bound on the rate of this decay. We will 
then clarify the difference in robustness between classi- 
cal chaotic systems and their quantum counterparts, in 
terms of the size of structures found in corresponding 
phase space functions (classical probability distributions 
and quantum Wigner functions). Finally, we will illus- 
trate the central issues with a numerical example, plac- 
ing bounds on the time needed for a (classically chaotic) 
quantum environment to decohere a quantum system of 
interest. 

Quantum lower bound for overlap decay and 
decoherence time. Using the projection operator 
P(t) — \ip2) (02 1 to rewrite the above-defined overlap as 
O q {t) = (ybi\P\%pi), and applying the Schrodinger equa- 
tion in the Heisenberg picture, we obtain 

where (■■•) = {ip\\ ■ ■ ■ Now, the uncertainty relation 
for V and P is AVAP > \([V,P])\/2, where (AV) 2 = 
(V 2 ) — (V) 2 is the variance of the operator V in the state 
0i, and similarly (AP) 2 = (P 2 ) - (P) 2 = O q -0\. Com- 
bining this with (0) gives 



dO q 
dt 



< 



d0 



dt 



< -AV(O q 



where the unitary evolution of 0jr (t) in the Hilbert space 
of the environment is generated by Hs+V^. The initially 



leading, after some algebra, to the inequality 
O q (t)> cos 2 AVdt^j =cos 2 0(i), 



(5) 



(6) 
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valid until <p(t) (which never decreases) reaches it/2. Note 
that by reversing the roles of ipi and 02 hi this ar- 
gument, we typically obtain a quantitatively different, 
though equally valid, result. We can therefore view AV 
appearing in (^|) as the spread of V in either state 0i or 
-02, whichever gives the tighter bound. 

When ipi and 02 represent states of a quantum envi- 
ronment (as discussed above) , then (||) gives the following 
lower bound on the decoherence times: 

T D > 7T?i/2AF, (7) 

where AV is to be interpreted as the typical value of AV^ 
during the decoherence process. 

Quantum and classical overlap in terms of phase 
space distributions. Apart from studying the sensitiv- 
ity of quantum evolution in its own right, we would like to 
compare it with classical sensitivity, particularly in the 
case of chaotic evolution. We will work with functions 
in phase space as these transparently suggest a classical 
counterpart of the quantum overlap O q (t). 

Equation |3| can be rewritten as |7j 

O q {t) = 2irh J W 1 (x,p,t)W 2 (x,p,t) dx dp, (8) 

where the W^s are Wigner functions corresponding to 
0i and ip2-, evolving under Hi and i?2- Let us now con- 
sider two classical phase space distributions, Li(x,p,t) 
and L2{x,p,t), obeying the Liouville equation under the 
respective classical Hamiltonians H\(x,p) and H2(x,p). 
Let us furthermore set the initial conditions for the Us to 
be the same as those for the W's: L\ = L% = W\ = W2 
at t = 0. In view of (||) it is now natural to define a 
classical overlap, 

O c (t) = 2nh J L 1 (x,p,t)L 2 (x,p,t) dx dp, (9) 

where the (arbitrary) normalization factor 2nh~ was cho- 
sen so that O c (0) = O q (0) = 1. By comparing the decay 
times of O q (t) and O c (t), we now have a setup for com- 
paring quantum and classical sensitivity to perturbations 
in the equations of motion. Because of the relevance to 
decoherence, we will refer to these as the quantum and 
classical decoherence times, though in the classical case 
this is just convenient nomenclature. 

The smallest structures of phase space distri- 
butions. A central hypothesis of this Letter is that the 
time scale for the decay of the overlaps O c and O q is de- 
termined primarily by the size of the smallest structures 
in the corresponding phase space distributions, with par- 
ticular relevance when the classical evolution is chaotic. 
In both cases we have two initially identical phase space 
distributions (the L's or W's) evolving with time while 
slowly accumulating a relative displacement due to the 
perturbation; a substantial decay of overlap occurs when 
this displacement is large enough that the two functions 



no longer "sit one on top of the other" . Clearly, this de- 
pends not only on the rate at which the functions move 
apart, but also on the local smallness of their structure, 
as this determines the degree of displacement needed to 
kill the overlap. The difference between the decay of 
overlap in the classical, chaotic case, and in its quantum 
counterpart, arises because of the qualitatively different 
mechanisms governing the emergence of small-scale de- 
tails in the corresponding phase space distribution. 

In the classical case, the size of local structure in 
L\ and L 2 shrinks exponentially with time, due to the 
stretching and folding associated with chaotic evolution: 
the probability distributions become thin and elongated, 
with a local width decreasing as exp(— Xt), where A is the 
largest Lyapunov exponent. As there is no lower bound 
to this smallness, it is clear that the decay time will be 
set predominantly by the Lyapunov time. By contrast, 
there are limits on the fineness of detail that can develop 
in the Wigner function, W; e.g. the Wigner function of 
a superposition of two identical Gaussians separated by 
AX - a Schrodinger cat-like state - exhibits interference 
fringes in momentum on a scale Sp ~ h/ AX 0. More 
generally, when spread over an area A — AXAP in two- 
dimensional phase space, W exhibits local structure on 
scales Sp ~ h/ AX, Sx ~ h/AP ||||. The corresponding 
phase space scale is associated with the sub-Planck ac- 
tion a ~ h 2 1 A which has physical consequences 0. Most 
notably (in the present context), the decay of O q occurs 
when the relative displacement of W\ and W2 is sufficient 
for their respective smallest-scale fringes to interfere de- 
structively. 

Two examples serve to build up intuition related to 
these issues. 

Example I. Let us assume that identical wave func- 
tions 01 and tp2 are superpositions of N Gaussians Gj ~ 
exp(— (x — Xj) 2 /(2h))e~xjp(ipjxh~) 

N 

0i (x) = 02 (x) = ^2Gj(x;xj, pj ) . 

3=1 

The corresponding Wigner functions W\ and W2 con- 
sist of N coherent-state Gaussians Gj , centered at points 
(xj,pj), as well as pairwise interference terms Gj^'. 

N 

3=1 j<k 

We assume (following ||) that the coherent-state Gaus- 
sians are sparse: each pair Gj and Gu is well separated 
by the "distance" dj^ phase space. The interference 
term Gj^k is then another Gaussian located halfway be- 
tween Gj and Gk, modulated by an oscillatory factor 
of frequency dj,k/Ti and twice the amplitude of Gj ||. 
The overlap (||) of W\ and W2 then works out to be 
O q « Ng + (N - l)Ng, where g = J G 2 dxdp. The first 
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term corresponds to the overlap deposited in the coherent 
Gaussians Gj, the second in the interference terms Gj^. 
Thus, for large N most of the overlap resides in the inter- 
ference terms. If we now displace one of the WF's relative 
to the other, by a distance at least on the order of the 
size of a typical interference fringe (but small compared 
to the size of the Gj's), then the contributions to the 
overlap from the Gj^s will typically interfere destruc- 
tively, resulting in a total overlap ~ 1/N <C 1. This is 
somewhat counterintuitive: for a fixed number of Gaus- 
sians of fixed size, we can increase the sensitivity of the 
system - as measured by the perturbation needed to kill 
the overlap (V'llV^)! 2 ~~ simply by increasing the average 
distance between the Gaussians, or equivalently the total 
area occupied in phase space. 

Example II. While classical probability distributions 
have no interference fringes, there is no bound on the 
smallness of structures resulting from chaotic stretching 
and folding. Let us examine a simplified model of the 
classical evolution of two Gaussian distributions in phase 
space, initially identical: 

1 _ * 2 _ p 2 

L 1 (x,p,t = 0)=L 2 (x,p,t = 0) = - e 

Zna x a p 

with x and p dimcnsionless. We now assume that with 
time both L's are exponentially stretched in the p di- 
rection and squeezed in the x direction, in an area- 
preserving way: cr p (t) = crexp(At), a x (t) = erexp(— At); 
furthermore, while the centroid of L\ remains fixed at 
x = p = 0, the centroid of L 2 drifts with a constant 
velocity v = (v x ,v p ). These assumptions mock up the 
relevant features of chaotic evolution under slightly dif- 
ferent Hamiltonians, where |v| is an indicator of the size 
of the perturbation. A simple calculation gives us the 
following decay of the overlap between L\ and L 2 : 

O c (t) = e -K*-e At /2.) 2 e -Kt-e--/2<x) 2 . ( 10) 

Generically, the first factor will dominate, and the over- 
lap will decay to negligible values on a time scale set by 
A" 1 . 

While both of these examples are highly simplified, 
we believe they capture the essential physics. We now 
present numerical results illustrating actual evolution. 

Numerical simulation. An example of a time- 
dependent Hamiltonian that generates chaos in one di- 
mension is [ fL0[ : 

2 ? 

p X 

H = Kcostx — Zsin(t)) + a— . (11) 

2m 2 

For parameter values m = 1, K = 0.36, I = 3.8 and 
a = 0.01, the stroboscopic Poincare surface of section, 
Fig.|l|, consists of four islands of stability surrounded by 
a chaotic sea. We have simulated both quantum and 
classical evolution, starting from a Gaussian distribution 



centered just outside the regular region (see Fig.|]) at 
t = 0, and evolving under H until time t = T, at which 
point a perturbation is turned on and the evolution forks 
into two branches governed by the Hamiltonians 

v 2 (x ± e) 2 

H± = - kcos(:e - Zsin(T + -r)) + a± '—, (12) 

2m 2 

where e = 0.5 and r = f — T. The perturbation is 
thus V = H + — H- = 2aex. The preparation time in- 
terval (0 < t < T) allows the distributions to develop 
small structures in phase space. After the perturbation 
is turned on at t = 0, we monitor the decay of overlaps. 




j I J : l, .. - - i □ I 

-10 10 

x 

FIG. 1. Poincare surface of section. The black ellipse shows 
the initial distribution. The evolution is confined to the dark 
cloud surrounding the islands of stability. 

Figure ^| shows how the decoherence time 7~d (defined 
here somewhat arbitrarily as the time r at which the 
overlap decreases to a value 0.9 Q) depends on prepa- 
ration time T for quantum and classical evolution. For 
short preparation times, both are equally sensitive to the 
perturbation applied at t = 0, reflecting the fact that 
the size of the smallest structure is basically the same in 
both cases. However, once the distributions have spread 
over much of the dynamically accessible area in the phase 
space, which occurs at T w 20, the size of interference 
fringes in the Wigner functions saturates, resulting in 
more or less constant decoherence times even for long 
T. Note that the quantum lower bound (^) denoted by 
crosses in the Fig.||, with AV evaluated directly from the 
simulation, gives results very close to the actual decoher- 
ence times; this indicates that the quantum states used 
in our evolution are close to minimum uncertainty states 
with respect to the uncertainty principle mentioned after 
(^). In the absence of such structure saturation in the 
classical distributions L, the classical decoherence time 
continues to decrease with increasing preparation time, 
due to the presence of ever smaller structures. While the 
computational cost of the classical simulations became 
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prohibitive for T > 30, the observed decrease in td is 
consistent with a rapid approach toward zero. 



classical 
Q—Q quantum 
X quantum lower bound 




FIG. 2. Time after which the overlap decreases to 0.9 ver- 
sus preparation time for classical and quantum environments. 

To further investigate the relevance of the smallest 
structure scale in the Wigner function to the decay of 
O q , let us define the spread of the system in position 
as AX = ((x 2 ) — (i) 2 ) 1 / 2 , and in momentum AP — 
((p 2 ) — (p) 2 ) 1 / 2 , with averages (■ • •) taken at t = 0. This 
translates to interference fringes of size Sx ss h/ AP in 
position and Sp m h/AX in momentum. Using our data, 
we have determined 5x and dp for each of the fourteen 
preparation times T shown in Fig.^|. 



that the value of 8p does not change much during the 
process. Hence, given a constant rate of relative drift in 
the momentum direction (due to the form of our pertur- 
bation, V = 2aex), we expect td oc Sp. On the other 
hand, when Sp is initially large (small T), then the deco- 
herence time will also be large, and 5p itself will decrease 
during this time; hence, we expect in this regime to ob- 
tain values of td that are smaller than suggested by the 
initial value of Sp, in agreement with the three highest 
data points shown in Fig.^. 

Conclusions. We have investigated the sensitivity of 
classical chaotic systems and their quantum counterparts 
to perturbations in their equations of motion. From gen- 
eral quantum considerations, we have derived a lower 
bound for the decay of O q (t). We have further argued 
that the sensitivity (in both the classical and quantum 
cases) is set by the size of the smallest structure of the 
related phase space functions. Finally, we have discussed 
the relevance of these results to the ability of quantum 
environments to rapidly decohere systems to which they 
are coupled. 

We are grateful to Diego Dalvit and Salman Habib 
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the Avalon Beowulf cluster at Los Alamos National Lab- 
oratory. This research is supported by the Department 
of Energy, under contract W-7405-ENG-36. 
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FIG. 3. Time after which overlap falls by 10% versus Sp. 

For the parameter values and form of the perturbation 
we have chosen, we have found that Sp is more relevant 
than Sx for the decay of O q (t). Therefore in Fig.||we plot 
the dependence of td on Sp. The linear dependence ob- 
served for small values of Sp (large T) is not unexpected. 
Recall that decoherence is achieved when the displace- 
ment of W2 relative to Wi becomes comparable to Sp. 
In the regime of small Sp, this occurs sufficiently rapidly 
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